A TRACE FORMULA FOR THE DIRAC OPERATOR 
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£*0 ' Abstract. Our goal is to extend the theory of the spectral shift function to 

the case where only the difference of some powers of the resolvents of self- 
adjoint operators belongs to the trace class. As an example, we consider a 
couple of Dirac operators. 
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Oh ■ 1. Introduction 

C/3 ■ 

The concept of the spectral shift function (SSF) first appeared in the work of I. M. 
Lifshits 12 in connection with the quantum theory of crystals. A mathematical 
theory of the SSF was shortly constructed by M. G. Krein in ^0]. One of his 
results can be formulated in the following way. Let Hq and H be self-adjoint 
operators with a trace class difference V = H — Hq . Then there exists the function 
£(A) = £(A; H, Hq), £ 6 Li(R), known as the spectral shift function (SSF) such that 
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the trace formula 

Tr(/(F) - f(H j) = I £(A)/'(A)dA, £(A) = £(A; H,H ), (1.1) 



^sO \ holds at least for all functions / s Cg°(R). Later in JT] M. G. Krein returned to 

this problem. In particular, he has shown that formula similar to ijl.lf) remains true 
for a couple of unitary operators Uq and U with a trace class difference. In terms 
of self-adjoint operators this means that formula (|l.l|l holds if the difference of the 
resolvents of the operators Hq and H belongs to the trace class Si (such operators 
Hq and H are called resolvent comparable). This allows to write formula (|l.lfl for 
sufficiently large class of differential operators. A relatively detailed presentation 
of the theory of the SSF can be found in |S] or ^S] . 

In 1957 [El 0E] T. Kato and M. Rosenblum have proven the existence of the 

S^ ■ wave operators (all relevant definitions can be found in J^l) for the pair of self- 

adjoint operators with a trace class difference. This result was extended in [2] by 
M. Sh. Birman and M. G. Krein to unitary operators. More important, in the same 
paper they have found a connection between the scattering theory and the theory 
of the SSF. Actually, they have shown that the corresponding scattering matrix 
S(X; H, Hq) minus the identity operator / belongs to the trace class and 

Det S(X; H, H ) = e r 2m ^ x ' H - Ha) (1.2) 

for almost all A from the core of the spectrum of the operator Hq. 

Later, in 8 T. Kato has proven the existence of the wave operators for the pair 
Hq, H under the assumption that 

R m {z)-Ro n {z) eSi, Imz^O. (1.3) 

However the construction of the SSF under this assumption seemed to be an open 
problem. Our goal here is to fill in this gap. Apart from its conceptual naturalness, 
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there are very simple applications which require such an extension of the theory. 
Here we give only one application - to the Dirac operator. 

We note that very general conditions of the existence of wave operators were 
obtained by M. Sh. Birman in the framework of the local trace class approach 
PP. On the other hand, the local theory of the SSF due to L.S. Koplienko |5] is 
somewhat less satisfactory. 

The author is grateful to M. Z. Solomyak for useful discussions of the Double 
Operator Integrals theory. 

2. The spectral shift function 
Let us first formulate the result of M. G. Krein obtained in jllj . 

Theorem 2.1. Let 

(h-t)- 1 -(ho-i)- 1 eSi. (2.1) 

Suppose that a function g(/J.) has two bounded derivatives and 

d a (g([i)-g n- 1 ) = O(\ri- 1 -^ a ), a = 0,1, 2, e > 0, 
where the constant go is the same for /i — > oo and fi — > — oo. Then 

g{H) - g(H ) G & ± 
and there exists the SSF £(//; h, ho) such that 

\ti(^h,ho)\(l + \n\r 2 dn<™ (2.2) 



(g(h) - g(h )) = f Z{n;h,h )9\lJ)dn. (2.3) 



Tr| 

Our goal is to extend this result to the case where condition 1(1. 3() is satisfied for 
some (not necessarily m = 1) odd m. 

Theorem 2.2. Let, for a pair of self-adjoint operators Ho and H , the assumption 
(11.3(1 hold for some odd m and all Imz ^ 0. Let a function /(A) have two bounded 
derivatives, and let 

d a (f(X)-foX- m ) = 0(\\\- m ^- a ), a = 0,1, 2, e>0, 

where the constant /o is the same for A — > oo and A — > — oo. Then the inclusion 

f(H) - f(Ho) G 6i (2.4) 

holds and there exists a function (the SSF) £(A;77, Ho) satisfying condition 

/oo 
|^(A; J ff, J ff )|(l + |A|)- m - 1 dA<oo (2.5) 

-oo 

such that the trace formula (11.11) is true. Moreover, for the corresponding scattering 
matrix S(X; H, Ho), the operator S(X; H, Ho) — I £ &i and relation (|1.2|) holds for 
almost all X from the core of the spectrum of the operator Hq ■ 

Our proof of Theorem 12 . 21 relies on its reduction to Theorem 12.11 although, sim- 
ilarly to |TI], we could have deduced it from the corresponding result for unitary 
operators. Actually, we construct a function ip such that the operators ho — (p(Ho) 
and h = ip(H) are resolvent comparable. To be more precise, we shall prove the 
following result. 



A TRACE FORMULA FOR THE DIRAC OPERATOR 3 

Theorem 2.3. Let the assumption 1(1. 3j) hold for some odd m and all Imz ^ 0. Let 

ip G C 2 (R) and let ip(X) — X m for sufficiently large \X\. Then the pair ho = ip(Ho), 
h = (p(H) satisfies condition (12. 1|) . 

We postpone the proof of Theorem 12.31 until the next section. Here we use it 
for the construction of the SSF. Since, under the assumptions of Theorem 12.21 
the operators ho = (fi(Ho) and h — (p(H) are resolvent comparable, we can apply 
Theorem 12.11 to the couple ho, h. Then we define the SSF for the pair Ho, H by 
the relation 

£(A; H, Ho) = S(<p(X);<p(H), <p(H )). (2.6) 

Suppose that the function ip is invertible and that 

ip'(X) > c> 0. 

Set fj, = >p(X), t/j = p^ 1 , g(/i) = f(ip(/J,)). Then formula ((2.3(1 implies that 

Tc(f{H) - f(H )) = Tr(g(h) - g(ho)) = I ^;h,h )g'(fi)dn. 



This coincides with l|l.l|) if the SSF £(A; H, H ) is defined by formula 

It follows from estimate 112.21 and the conditions on <p(X) that the function l|2.6|) 
satisfies estimate (|2.5|) . The class of functions /(A) for which formula Hl.lfl is true 
is obtained from the class of functions g(fi) by the change of variables fi = y(A). 
Finally, formula (|1.2|l follows from the same formula for the pair ho, h, definition 
(12.61) of the SSF and the invariance principle for scattering matrices. This concludes 
the proof of Theorem 12 . 21 given Theorem 12.31 

We emphasize that the trace formula l|l.l|) fixes the SSF £(A; H, Ho) up to an ad- 
ditive constant only. This constant remains undetermined also by condition l|2.5fl . 
Nevertheless (see ^J]) defining the SSF via the corresponding perturbation deter- 
minant fixes £(X; H, Ho) up to an integer constant. It is exactly for this "correct" 
choice of a constant that relation 11.2(1 holds. 

3. Proof of Theorem 12.31 

Our proof of Theorem 12.31 relies on the theory of Double Integral Operators 
(DOI). Let us briefly recall its basic notions (see |3II1]; for details). Let Ho and H 
be a pair of self-adjoint operators in a Hilbert space 7i. Denote by Eq and E their 
spectral families. Let us use the following result which can be deduced from 0] . 

Proposition 3.1. Let 

/oo />oo 
/ K(X,fi)dE(fi)TdE Q (X), (3.1) 

-oo J — oo 

where the kernel K(X,n) is bounded, i.e., 

\K(X,fi)\<C <oo, (3.2) 

it is differentiable in X and 

\dK(X,fi)/dX\<C{l + X 2 )-\ (3.3) 

Assume, moreover, that 

lim K(X,n)= lim K{X,p) (3.4) 

A— » + oo A— > — oo 
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(these limits exist by virtue of Q3.3JI ). Then the transformer $ : Si — > 61 defined 
by H3.1|) is bounded. 

Let us set 

5 (A)=g 2 (A) = (A-z)- m 
and 

T = T z = g z {H)-g z (H Q ). 

We need the representation of the difference /(if) — f(Ho) hi terms of the DOf : 

/(#) - /(jy ) = /°° r M^ffi( P )M„(A). (3.5) 

J-ooJ-oo 3(A) -3(W 

A technical problem with a proof of Theorem l2 . 3l is that the denominator g(X)—g(fi) 
in (|3.5|) has "extra" zeros on the antidiagonal A = — /1. Therefore we shall split 
Theorem 12.31 into two separate assertions. 

Proposition 3.2. Let f G Cq(R), and let supp/ C [— r, r\. Let the assumption 
(|1.3I) be satisfied for some odd m and z = ia with a sufficiently large (compared to 
r) a. Then inclusion (|2.4() holds. 

Proposition 3.3. Let E C 2 (R), 0(A) = for |A| < r and 0(A) = 1 for |A| > 2r. 
Set 

/(A) = 0(A)(A m -z)- 1 . (3.6) 

Then inclusion Ij2.4|l ZioWs i/ the assumption H1.3(l zs satisfied for some odd m and 
z = ia with a sufficiently small (compared to r) a. 

Propositions 13. 21 and 13. 31 imply of course Theorem l2.3l Indeed, 

(<p(H) - i)- 1 - (<p(H ) - i)- 1 = (f Q (H) - /„(#„)) + (f(H) - f(H )), (3.7) 

where /o(A) = (1 — 9(A))((p(\) — i) _1 has finite support and / is given by formula 
H3.bjl . Both terms in the right-hand side of (|3.7|) belong to the trace class. 

For the proof of Propositions l3 . 21 and !3 . 31 we need the following elementary result. 

Lemma 3.4. Set 

p(\, w z) = (A - z)" 1 " 1 + (A - z) m " 2 (/i - z ) + ... + ( Ji - z ) m -\ (3.8) 

Let m be odd, and let r be some fixed number. Then for |A| < r, |/i| < r, z = ia 
and a sufficiently large a 

\p(X,n;z)\>c>0. (3.9) 

Similarly, if either |A| > r or \fi\ > r, z — ia and a is sufficiently small, then 

\p(\^;z)\>c(\\\ + \fi\) m -\ c>0. (3.10) 

Proof. By virtue of the equality 

p(\, (i; z) = (X- ia)™- 1 ^ + 0- + --- + cr™" 1 ), o = (ji- ia)(\ - ia)' 1 , 

both estimates l|3.9|l and (|3.1U|) (by the proof of H3.1(J[) we assume that \fi\ < |A|) 
reduce to the same estimate 

|1 + (T + ••• + cr m - 1 \ > c>0. (3.11) 

If |A| < r, \fi\ < r and a is sufficiently large, then 

a = (1 + i/ia _1 )(l + iAa -1 ) -1 

belongs to a neighbourhood of the point 1 which implies l|3.11|) . 
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Remark that the zeros of the function 1 + a + ■ ■ ■ + a" 1 ^ 1 are given by the formula 
(7k = cxp(27tik/m), k = 1, . . . , m — 1. 
If |A| > r, |A| > |/i| and a is sufficiently small, then the values of 

a = (x — ie)(l — is)^ 1 , x = /iA _1 e [—1, 1], e = aX^ 1 , 

belong to a neighbourhood of the real axis and therefore are separated from all 
points (Tfc. This again implies (|3.11l) . □ 

For the proofs of Propositions 13 . 21 and 13.31 we consider DOI l|3.5|l with kernel 

m,) = ^-44 (3.12) 

and verify the assumptions of Propositions 13.11 Obviously, the function K(X, /i) 
tends to f{n)g{^)~ 1 as A — > ±oo so that condition l|3.4|) is satished. The estimates 
(13.21) . 13.3fl we shall check separately under the assumptions of Propositions l3 . 21 and 
13.31 Note previously that, under the assumptions of Proposition l3.2l K(X, /i) = if 
|A| > r and \^i\ > r. Similarly, under the assumptions of Proposition l3~!fl K(X,fi) = 
if |A| < r and |/i| < r. 
Since 

ff (A) - g( P ) = (A - z)- m (p - z)- m (n - X)p(X, n; z), (3.13) 

we have that 

K(X,n) = -$(X,u)G(X,n), (3.14) 

where 

*{\,ri= m - f{fi) (3.15) 

A — /! 

and 

(A-zH/i-z)" 

G(A,m) = tt r . 3.16) 

p(X,n;z) 

Let us start with Proposition 13.21 and consider first the region where |A| < R, 

|/i| < R for some R. Here we proceed from representation (|3.14(l . The function 

(|3.15J) is bounded because f £ C 1 . Since 

|/(M) - /(A) - /'(A)(/i - A) | < 2- x (/x - A) 2 supp M < fi |/»|, 

the function 

/(M)-/ (A)-f( A)( M -A) 

(A-M) 2 

is also bounded for |A| < R, |/i| < R. According to Lemma 13.41 the function (|3.1t)|) 
as well as its derivatives in A are bounded in this region provided z = ia and a is 
sufficiently large. 

It remains to consider the region |/i| < r, \X\ > R > r (and, similarly, |A| < r, 
\fi\ > R > r). Here we use that the denominator in (|3.12|> is bounded from below 
because 

\gW-9(ji)\ > {p 2 +a 2 )- m/2 -{X 2 +a 2 r m/2 > (r 2 +a 2 )- ,n/2 -(R 2 +a 2 )- m/2 > c> 0. 

(3.18) 
Therefore the function (|3.12|) is bounded. Further, differentiating l|3.12|l . we find 
that 

dK(\n) _ /'(A) (/(A) - /MMA) 



9$(A, „)/dX = MW JV ;i =Lf^ " ; . (3.17) 



5A g(\)-g(ji) (g(X)-g(v)y 



(3.19) 
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Thus, it follows from (|3~TS|) for |/x| < r, |A| > i? > r or for |A| < r, |/x| > i? > r that 

^'^ l < Cd/'WI + (l/WI + l/(M)l)l</(A)|) < d(l + IAD— 1 . (3.20) 
This concludes the proof of Proposition 13. 21 

Now we check the estimates lKS.2t and IJ3.3I1 under the assumptions of Proposi- 
tion 13.31 

Let us first consider the region where |A| < 3r, \/x\ < 3r (with the square |A| < r, 
|/i| < r removed). Similarly to the proof of Proposition 13.21 we have that func- 
tions l|3.15(l and (|3.17(l are bounded. The boundedness in this region of the func- 
tion i|3.16[l . as well as of its derivatives in A, follows, for sufficiently small a, from 
Lemma 13.41 

Next we consider the region |A| > 2r, \fj,\ > 2r where 

/(A) - /( M ) = (X m - i)- 1 - ( M m - z)- 1 = (\ m - i)-\n m - _1 (a* - a)p(A, M ; o) 

and the function p(A, fx; 0) is defined by formula (|3.8(l . Hence it follows from (|3.12(l 
and i(3~T3|> that 

K(\ U ) = (A-^r(M-^XA,/i;0) 
1 ' X m -i n m -i p(\,n;z)' 
By virtue of estimate (|3.10|l . this function is bounded. For the proof of (|3.3() . we 
use two estimates 



d {X-z) 



dX X m 
and 

0*»(A, W O) 



<cx- 2 



dXp(X, n; z 
The first of them is obvious. For the proof of the second, we remark that 

d p(X,n;0) 



<C(|A| + H)- 2 . (3.21) 

mark tha 



dXp(X, ix; z) 

((p(A, M5 z ) ~ P{\ M5 0))p' x (X, ft z) - (p' A (A, fi; z) - p A (A, ^; 0))p(A, m; z )J • (3-22) 

Let us now take into account that p(A, /i; z) is a polynomial in A and /i of degree 
m — 1, but the terms of order to — 1 in p(X, /i; z) and p(X, \x\ 0) cancel each other. 
Therefore their difference consists of terms X p fi q where p + q < m — 2. This gives 
the estimate 

|p(A, /i; z) - p(X, n; 0)| < C(|A| + H)" 1 " 2 . (3.23) 

Similarly, differentiating the difference p(X, /i; z) — p(X, n; 0), we see that 

d(p(X, m; z) - p(X, m; 0))/9AJ < C(|A| + H) m - 3 . (3.24) 

Substituting (|3~2l|l and (J3~2lj> into (JS23, we obtain lt3~2lTl . 

Let us finally consider the region where |A| > 3r, |/x| < 2r or |/z| > 3r, |A| < 2r. 
Here we use again representation (|3.12(l and estimate (|3.18|l which yields 

|tf(A,/0|<C(|/(A)| + |/(/K)|)<Ci<oo. 
The same bound (J3.18JI shows that the derivative H3.19|l satisfies estimate H3.2()|l . 
which yields l|3.3(l . 

This concludes the proof of (|3.3(l and hence that of Proposition EH 
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4. The Dirac operator 

As an example to which Theorem 12 . 21 directly applies, we now consider the Dirac 
operator describing a relativistic particle of spin 1/2. Let Ti = L2(M. 3 ; C 4 ) and 

3 

Hoo = ^2 a i D i + ma °' ( 4 ' 1 ) 

3=1 

where m > is the mass of a particle and 4x4- Dirac matrices satisfy the 
anticommutation relations 

atiOtj + otjQti = 0, i 7^ j, oq = I. 

These relations determine the matrices a, up to a unitary equivalence in the space 
C 4 . Their concrete choice is of no importance. If Qij are replaced by uatiu* where 
u is a unitary transformation in C 4 , then the operator Hoo is replaced by a unitary 
equivalent operator of the same structure. 

Making the Fourier transform <£>, we find that Hoo — <&*A<i> where A is multipli- 
cation by the matrix function (the symbol of iJoo) 

3 

A(g) = ^atjZj +ma . 

3=1 

It is easy to see that A(£) has the eigenvalues 

a h 2(0 = -a3A(0 = (\e + ™ 2 ) 1/2 (4-2) 

of multiplicity 2 so that 

A(0 = T(0A(0T*(0 (4.3) 

where the matrices T(^) are unitary and A(£) = diag{oi(^), 02(^)1 a 3(0i°4(C)}- ^ 
concrete form of the matrices T(£) is inessential for us. In particular, the operator 
Hqo is selfadjoint on the Sobolev space H (K 3 ; C 4 ) =: V(Hoo). 
Here we consider a couple of Dirac operators 

H = Hoo + V , H = H + V (4.4) 

where Vo is multiplication by a symmetric bounded 4x4- matrix function Vq(x) 
and a perturbation V^ which is also a symmetric 4x4- matrix function satisfies 
the condition 

\V(x)\<C(l + \x\)- p , P >3. (4.5) 

We do not make any special assumptions on matrices Vo(x) and V(x). In particular, 
the spectrum of the operator H might cover the whole real axis. Particular cases 
3 3 

V (x) = ^2vf\x) aj +v {0 \x), V(x) =^v {x) aj +v{x), 

3=0 J=0 

where in \v^ and Vj , v are scalar functions, correspond to an interaction of an elec- 
tron (or of a positron) with magnetic and electric fields (with potentials (v%, v%, V3) 
and v, respectively). We suppose that the "background" potential Vq(x) is only a 
bounded function whereas the perturbation V(x) satisfies the condition 14.5fl . We 
shall show that in this case condition (|1.3|l is satisfied for m = 3. 
Differentiating the resolvent identity 

R(z)-Rq(z) = -R(z)VRq(z), Imz ^ 0, 
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we find that 

m 

R m {z) - R™(z) = - ^ R k {z)VR™ +1 - k {z). (4.6) 

fe=i 
If we replace the resolvents R and i?o by i?oo in the right-hand side of (|4.6|l , then 
we obtain the sum of terms Rq VR™ +1 ^\ It is easy to see (this follows from 
Proposition 14. II below ) that, for m > 3, these operators belong to the trace class. 
Therefore it suffices to justify the replacement of R and Rq by i?oo- However the 
boundedness of the operator (H 00 — z) k R k {z) is equivalent to the inclusion 

V(H k ) C V(H k ). 

If k > 1, this inclusion requires boundedness of derivatives of the function V and is, 
in general, violated under the conditions above. To bypass this difficulty, we suggest 
a trick based on commutation of the operators (x)~ r and R(z). We introduce also 
the whole scale of symmetrically normcd ideals & p . 
We start however with operators (x)~ t Rq (z). 

Proposition 4.1. // 

p > 3/min{r, k} =: p(r, k), p > 1, (4-7) 

then 

(x)- r RUz) G © p . 

Proof. Since i?oo( z ) = < &*(^4 — z)~ fc $, we have to check that the integral operator 
with kernel 

(xr r eM*<x,t>)(MO-z)- k 
belongs to the class 6 P with p determined by l|4. 7fl . According to l|4.2|l . (|4.3|) . the 
function 

(e) fe (A(o-z)- fc 

is bounded so that it suffices to consider the operator with kernel 

(a;)- r exp(z<x,C>)(e) _fe . 

This operator belongs to the required class & p according, for example, to the results 
ofJEj. □ 

Next we extend this result to the operator H (or Hq). 

Proposition 4.2. If p satisfies Ij4.7|l . then 

(x)- r R k (z) e 6 P . 

Proof If k = 1, then 

(x)- r R(z) = ((x)- r Ro (z)) ■ {{Hqo - z)R(z)) 

and the second factor in the right-hand side is a bounded operator. Let us justify 
the passage from k to k + 1. Since 

(H - z)(x)- r - (x)' r (H -z) = [tf 00 , (x)- r ], 

we have that 

(x)- r R k+1 = R(x)- r R k + R[H n, (x)- r }R k+1 . (4.8) 

Let us write the first term in the right-hand side of (|4.8|l as 

R(x)-' r R k = {R(x)- r «){(x)- kr «R k ) (4.9) 
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where r = r(k + l) -1 . Here R(x)~ r " e & p for p > p(r Q , 1), and (x)~ kr °R k £ & p 
for p > p(kro, k). Therefore the product (|4.9(l belongs to the class & p where 

p- 1 < p(r , l)' 1 +p{kr , k)- 1 = p{r, k + l)" 1 . 

The second term in the right-hand side of (|4.8|1 is even better since the matrix 
[-ffooi ( x )~ r ] is bounded by (x)~ r ~ 1 . D 

Now it easy to verify inclusion 1(1.3(1 for the Dirac operators. 

Theorem 4.3. Suppose that the symmetric matrix functions Vq(x) and V{x) are 
bounded and V{x) satisfies condition (|4.5(l with p > 3. Then condition ((1.3(1 is 
satisfied for all m > 3. 

Proof. Indeed, let us write all terms in the right-hand side of (|4.(i|) as 

R k VR ™+i-k = (R''(x)- r i)((x)' > V)({x)- r *RZ +1 - k ) (4.10) 

where ri = kp{m+ 1) , 7"2 = (m + 1 — k)p(m + 1) . According to Proposition ^. 21 
the first term in the right-hand side of 1(4.10(1 belongs to the class 6 P for p > p\ = 
p(ri, fc) and the last term belongs to the class & p for p > p2 — p(r2,m + 1 — k). 
Since p^ 1 + p^ 1 > 1, the product l|4.10|l is trace class. □ 

Corollary 4.4. The WO W±(H,Hq) exist and are complete. 

Finally, combining Theorem 14 . 31 with Theorem 12. 21 we obtain 

Theorem 4.5. Let the operators Hq and H be given by formula Ij4.4|l where Hqq 
is the "free " Dirac operator l|4.1|l . Suppose that the symmetric matrix functions 
Vq(x) and V(x) are bounded and V(x) satisfies condition (|4.5|l with p > 3. Let a 
function /(A) have two bounded derivatives and 

d a (f(\)-f \- 3 ) = O(\\\- 3 - e - a ), a = 0,1, 2, e>0, 

where the constant /o is the same for A — > oo and A — > — oo. Then inclusion 
(12.41) holds and there exists a function (the SSF) £(A; H, Hq) satisfying condition 
(12.5(1 where m — 3 such that the trace formula (jl.ljl is true. Moreover, for the 
corresponding scattering matrix «S*(A; H, Hq), the operator 5*(A; H, Hq) — I € @i anrf 
relation ((1.2(1 ZioWs /or almost all A /rom f/ie core o/ the spectrum of the operator 
Hq. 
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